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Abstract 

Using charged hard spheres model as an example, the dense one-component 
plasma is considered. For this model the Enskog-Landau kinetic equation is ob- 
tained and its normal solution is found using Chapman-Enskog method. Transport 
coefficients are obtained numerically and analytically and compared with the ex- 
perimental data available. 
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Construction of kinetic equations for dense gases and plasma is one of the most important 
problem in the kinetic theory of classical systems. A consequent approach for construction 
of kinetic equations has been done by Bogolubov [0. This approach is based on a chain of 
equations for s-particle distribution functions and on boundary conditions for weakening 
correlations. Using such boundary conditions, we can in principle, express all s-particle 
distribution functions in terms of the single-particle function and obtain for it a closed 
kinetic equation. There is a large number of approaches for derivation of kinetic equations 
|| ||, f|, H|. Despite a difference of these approaches between themselves in shape, the 
weakening correlation principle in one or another form has been used in all approaches 
just as in Bogolubov's method. However, all these approaches are most efficient in the 
case when a small parameter (density, interaction, etc.) is present. For dense gases and 
dense plasma small parameters are absent. In this case analysis of the BBGKY hierarchy 
becomes very difficult because we can not restrict ourselves to some finite number of 
terms in expansion for the collision integral. Moreover, an additional problem, concerning 
correct account in the collision integrals of interactions between particles on short as well 
as long distances, arises. 

Relatively recently an approach, based on a modification of the weakening correlation 
principle, has been proposed |J. This approach leads to a consequent construction of 
kinetic equations for dense gases without additional phenomenological assumptions. New 
boundary condition to the BBGKY hierarchy take into account a non-equilibriumnes of 
single particle distribution function as well as local conservation laws of mass, momentum 
and energy, i.e., the quantities which constitute the basic for the hydrodynamic description 
of evolution of the system. In the "pair collision" approximation, such approach leads 
to an Enskog-like kinetic equation. Similar ideas have been proposed independently by 
Karkheck, van Beijeren, de Schepper and Stell at derivation of the kinetic equation for 
the "square-well" potential. Somewhat different modification of Bogolubov's approach 
has been considered by Rudyak || ||. Here, the Enskog-like kinetic equation for a system 
of hard spheres has been obtained and attempts to extend this equation on a system with 
soft potential have been made. 

The ideas of work ||, which is based on Zubarev's non-equilibrium statistical operator 
method |L0], |ll|] , stimulate a revision of the problem connected with constructing of kinetic 
equations for dense gases and plasma. Investigations J12, |H| were logical continuation of 
the work || in which original result has been achieved: a consequent derivation of the 
kinetic equation of revised Enskog theory |14|, 15, pij] for a system of hard spheres. 

In the present paper a kinetic equation for the single-particle distribution function 
is obtained from the BBGKY hierarchy with modified boundary condition in the "pair 
collision" approximation. This kinetic equation is valid for moderately dense classical 
systems with the interparticle potential in a form of hard sphere potential plus some 
long-range potential <P'(r). In the case when <P l (r) is the Coulomb potential, we have 
obtained a kinetic equation, called Enskog-Landau one, for a system of charged hard 
spheres. Normal solutions of this equation are found by the Chapman-Enskog method. 
On the basis of the solutions, analytical expressions for viscosity coefficients and thermal 
conductivity are obtained. Numerical calculations of transport coefficients are performed 
for neutral and ionized argon. The results are presented in a temperature dependent 
form. A comparison between theoretically predicted values for transport coefficient and 
experimental data is examined. 
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The BBGKY hierarchy of equations for non-equilibrium distribution functions of classical 
interacting particles has been obtained in the paper || on the basis of assembling time 
retarded solutions for Liouville equation with modified Bogolubov's condition meaning 
weakening correlations between particles. According to Zubarev's non-equilibrium statis- 
tical operator method JlO, [11], full non-equilibrium distribution function p (x N ;t) for all 
N particles of the system satisfies the following asymptotic condition: 

t Jim^exp (tL N t ) (p (x N ,to) - p q (x N ,t J^ = 0. (1) 

Here, the limit to — > — oo is made after thermodynamical one N — > oo, V — > oo, 
N/V — > const, i = and is the Liouville operator: 



N I N N 

3=1 Z j=l k=l 



2m dr^ ' drj k \dpj dp k J ' 

<&jk is the interaction energy between two particles j and k; Xj = {r, p] is the set of phase 
variables (coordinates and momenta). Quasi-equilibrium distribution function p q (x N ;tj 
is determined from the condition of maximum for informational entropy at fixed values 
of the single particle distribution function fi {x x \ t) and average density of the interaction 

energy (£ int («*))*, . .)' = J dT^ . . . X p (x N ; tj, dT^ = - ^ j , that corresponds to 

taking into account correlations, related to conservations laws of hydrodynamical variables 
for particle density n (r; t), momentum j (r; t) and full energy £ (r; t) |L7j]. This function 
can be presented as follows || |13| : 

p q (x N ,t) = exp (-U N (r N ;t)) f[ (3) 

where u (r 3 -; £) is obtained from the relations: 



u(ri,t) 



J (JV-1)! v v JJ f± 2 u{rj] 



ty 



N j 

^tv t) = E ^ifeft*) ft* = r ^ *) = o *) + £ *) ) ' 

n(r;t) = J dpfi(x;t) is non-equilibrium particles concentration, /3 is certain function, 

being an analogue of local inverse temperature. 

Taking into account the boundary condition (HD is equivalent to the transition from 



the Liouville equation to a modified one [11 



+ ^ p (x N , t) = -e(p (x N , t) - p q (x N , t) ) . (4) 

This equation contains the small source in the right-hand side, which destroys the invari- 
ance with respect to time inversion (e — > +0 after the thermodynamic limit transition). 
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chain for the s-particle non-equilibrium distribution function f s (x s ; t) — J dT N _ s p (x N ; 
i, HI: 



(/>(•''■:/) .</.- (r s ;t)n/i(^;t)|, 



(5) 



where 



g s (r s ; t) = J dT N _ s dp s p q (x N , t) (6) 

is the quasi-equilibrium s-particle coordinate distribution function which depends on 
n(r;t) and f3(r;t) functionally. Due to the fact, that g± (r^t) = 1, the equation chain 
(j5|) is distinguished from the ordinary BBGKY hierarchy Q by the availability of sources 
in the right-hand parts of the equations beginning from the second one and takes into 
account both the one-particle and collective hydrodynamical effects. 

Let us consider the solution to equation chain (|j) within the pair collision approxima- 
tion. In this case, three- and higher-particle correlations are neglected, but environment 
influence to the evolution of a chosen pair of particles is taken into consideration by means 
of application of the boundary condition. Then from (|) for f 2 (xi,X2]t) we obtain |J: 

+ iL 2 + ej f 2 (x 1} x 2 ; t) = eg 2 (r u r 2 ; t) f\ (x X ] t) f l (x 2 ; t) . (7) 

The formal solution of equation (0) is of the form 

h (xi,x 2 ;t) = 

o 

= e J dTexp((e + tL 2 )T^g 2 (r 1 ,r 2 ;t + T)f 1 (x 1 ;t + T)f 1 (x 2 ;t + T). {b) 

— oo 

Substituting (|) into (H) at s = 1 leads to the kinetic equation for one-particle distribution 
function in the pair collision approximation: 

\^+tL(l)jf 1 (x 1 ;t)=I ca (x 1 -,t), (9) 

where 

o 

Icoi(xi,t) = - J dx 2 tL(l,2)e J drexp ( (e + iL 2 ) t)F 2 (xi, t + r) , (10) 

— oo 

F 2 (xi, t + t) = g 2 (n, r 2 ; t + r) f l (x x \ t + r) f x (x 2 ; t + r) 

is the collision integral. We must emphasize that equation @ is needed to be adjusted 
with the equation for binary quasi-equilibrium distribution function g 2 {r 1 ,r 2 ;i). Ac- 
cording to (0) and (|6|) this function functionally depends on ?i £i n t(y\t) (or on 
(3(r;t)). Besides it was shown, that quasi-equilibrium correlation distribution function 
h 2 r 2 ; t), which is related to g 2 (ri, r 2 \ t) (h 2 = g 2 — 1), satisfies the Ornstein-Zernike 
equation |12| . In paper [18] a non-equilibrium grand canonical distribution for the system 
of hard spheres and on the base of methods of non-equilibrium statistical mechanics [l9j an 
Ornstein-Zernike equation for pair quasi-equilibrium correlation function of hard spheres 
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paper |12| . New equation for h 2 (ri, r 2 ; t) is an analogue of Ornstein-Zernike equation 
in equilibrium statistical mechanics |19[ . This equation for the hard sphere system in 
equilibrium case has exact solution in Percus-Yevick approximation [0 . 

In the paper |6|] some special cases were considered, when the interparticle interaction 
potential <Pkj is modelled as the hard sphere potential <P hs (|r|) for particles with diameter 
a. Then taking into account the singularity of the hard sphere potential (r — > + 0, r 
is collision time) it was shown for the first time how the collision integral (|10D transforms 
into collision integral of revised Enskog theory (RET) 



1% (*i; *) = / dx 2 f hs (l, 2)g* s (r ls r 2 ; t) }\ (x l5 1) f x (x 2 ; t) 



where T hs (l,2) is the Enskog's collision operator for hard spheres, g 2 s (t*i, r 2 ; £) is the 
pair quasiequilibrium distribution function of hard spheres, which depends on the average 
non-equilibrium density n (r; t) functionally. It is important to emphasize that if -theorem 
for the kinetic equation (|]) with the Enskog collisions integral (|TT| ) has been proved by 



P.Resibois 15, 16 



In the same paper || another case was also considered, when the interparticle interac- 
tion potential is modelled as a sum of a short-range potential (hard spheres, for example) 
and some long-range smooth potential; 



where a* is the effective diameter of hard spheres, which depends on the method of 
splitting the potential $(H) into short- and long-range parts. 

If the time retarding and spatial inhomogeneity will be neglected, then we can present 
collision integral fllTf ) in the second approximation with respect to interacting potential 
<P(\r\) as follows pfl: 



J 0fe 


\r\ < a* 


1 


\r\ > a* 



I col {x x ; t) = (xi; t) + Ci fa t) + I l co i {xil *) , 
C (^15 t)=J dx 2 f hs (l, 2)g 2 (r 1; r 2 ; t) h fa; t) fi (x 2 ; t) 
iZi fail *) = / dx 2 tL l (l, 2)g 2 (n, r 2 ; t) / 2 (x x ; t) /i (z 2 ; t) 



4/ : 
1 d 



(12) 
(13) 
(14) 

(15) 



m dvi , 



d 



d 



dv is dv 2 „ 



fi {xi;t)fi (ri+r 12 ,v 2 ;t) 



where 



(16) 



1 °° 

— / dr 12 / da r\ 2 g 2 (n,n + r 12 ;t) 
m J J 



and 



g = V 2 - Vi, 



~d<P l (\r 12 \Y 


jdr 


~d<P l 


{\ri2+gr\)' 


dr 12)a 


— CO 




dr 12tf3 


= (r 2 - n 


)/\V2 ~ 


r i\- 





The first term in the right hand part of ([12]) is the Enskog collision integral (p~3|) , where 
92 r 2 ; t) is the quasi-equilibrium pair distribution function for system of particles with 
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and inverse temperature (3 (r; t) functionally. The second term (|14D is the mean field 
influence, and the third one ( |l"5|) is written in the form of Landau-like collision integral. 

If ^(|t*|) = - — — is the Coulomb potential, equation ( |l"2"D with structure of (|T6|) is a 
generalization of Landau collision integral without divergency for short-range distances 



between particles (that is different from the usual Landau collision integral [y, ^0] because 
short-range interactions are taken into account correctly). Therefore we can call such 
kinetic equation (|9]) with collision integral (0) as the Enskog-Landau kinetic equation for 
the system of charged hard spheres. It is necessary to emphasize that in equation (|l6|) 
the long-range divergency still remains. To avoid this problem sequentially we have to 
consider kinetic equation with taking into account of dynamical screening effects |2|, [3|. 
But this way is impossible in Enskog-Landau kinetic equation. Only one we can do for 
further calculation is to change upper integral limit to some finite value, which could have 
a meaning of value of statical screening in our system (see below), solve this problem we 
must consider dynamical screening effects. 

Following JTJJ, we can write for (xi;t) the next representation: 

C (*i; t) = IT (x 1 ; t) + I h JV ( Xl ; t) , (17) 

-C 0) ( x ^ f ) = / dv 2 debdbgg 2 [a + \n(r] t), /3(r; t)) x 



x [fi {ri, v[; t) fi (n, v' 2 ; t) - }\ (n, vt] t) f\ (r h v 2 ; t) 
^f 1) ( a; i;0 =cj3 / drisdv 2 (r 12 g) © (r 12 g) r 12 g 2 (cr + |n(r;t), (3{r;t)) x 

x (/i (r Xj «i; t) V/i (r ls v' 2 - 1) + h (n, w i; t) V/i (r ls v 2 ; t) ). (19) 

Here we use definitions: e is an azimuthal angle of scattering, f 12 is the unit vector, 
g — v 2 — Vi, b is the impact parameter, 0(x) is the unit function, i/j = v\ + f i 2 (f 12 • g), 
v 2 = v 2 — fi 2 (f 12 ■ fl 1 ) are velocities of particles after a collision, and cr + is a diameter of 
the particle. 

By representing (|l"5|) in the cylindrical coordinates, one can introduce the impact 
parameter b, azimuthal angle of scattering e, distance along the cylinder axis £. Then 
Boltzmann-like collision integral can be obtained from (p~5|) , putting g 2 = 1: 



J L *) = / ^2 cfe & dbg(fx (r*i, u*; £) /i (n, u 2 ; t) - /i (r x , «u £) A (n, v 2 ; t) ), (20) 



where 

= «i + Au, v% = v 2 — Av, 



i + ~ ( 21 ) 

Au = / ^V$ z (|r 12 |) 

ma J 



After these transformations we have the kinetic equation (|9|), where the collision inte- 
gral I co i(xi,t) can be presented as the sum of (|14j) , (]T8|) , ([19]) and (p0|). 
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We shall solve the equation (^) by means of iterations. Therefore a set of hydrodynamical 
variables should be introduced for reduced description of the system: hydrodynamical 
density of mass, momentum and kinetic energy j|, |2l| . The conservation laws are to be 
written down after multiplication of both left and right parts of equation @ by com- 



ponents of additive invariant vector m, mv as well as by kinetic energy 



mc 



and after 



integrating over v |4], [21| . The local-equilibrium Maxwell distribution function may be 
taken as an initial approximation: 



A (0) (xi;t) = n(n;t) 



in 



3/2 



2nkT (n;t) 



exp 



mc\ (t*i; t) 
'2kT(r 1 ;t) i 



(22) 



where n(r;t) is the density, c(r*;i) = v — V(r;t), and V(r;t) is the hydrodynamical 
velocity. The total distribution function fi(x;t) has a form /°(a;;i)fl + (p(xi,t)^j and 
the correction y?(xi;t) may be expressed through Sonine-Laguerre polynomials The 
calculations show lOl: 



Ail) 



m 



/5 mcf (ri; t) ' 
T(r i; t)V 2kT(n;t) \2 ~ 2fcT(n;t) / 

m5 (0) / 1 



(ci.V)T(n;*)- 



c 1 c 1 --4 I) : (Vy(n;t)) 



(23) 



2kT(ri,t) V x " 3 

where / is the unit tensor; A(l), -B(O) are coefficients which satisfy the following relations: 

1 + fvmcr 3 ^ (<7 + |n,/3) 



Al 



15 /7T 

¥V2 



X 



b(0) = 5 



7rm 



x 



n(a, (a+K/3)fi£ 2) + fif' 2) 

1 + ^nno z g 2 (a + \n 1 j3) 
n(^(a+|n,/?)fi£ 2) +n I (2 ' 2) ) 



(24) 
(25) 



Here, we use notations: 







(p.?) 

hs,l 



dgogl q+3 exp (-# 2 ) fig 

£7 

fig = 2n J bdb(l - cos V 





oo 

2tt / 1 



cos P x* 



00 



2fcT 



9 



(26) 

(27) 
(28) 



The expressions ( pq - |2q ) are known as fi-integrals @], x'i X* being the angles of scattering 



for the hard spheres and Coulomb particles respectively. 

The fi-integrals can be calculated exactly @] or approximately, if the first way fails. 
We shall attempt to do this from the geometrical point of view. Considering the dynamics 
of collision of hard spheres, we have [|| 



cos ■ 



X 



b 

2' 



and 
7 



fi 



(2,2) 
hs 



2vra 2 . 



(29) 



imation can be obtained [[22 



sinx" ~ 2'—^- I — — r - (30) 

m 9o J (b 2 + £ 2 )i 

(2 2) 

This expression for sinx* leads to logarithmical divergency in Q\ ' at integration over 
infinite sight parameter (|28|) . To avoid this difficulty, the infinite sight parameter in the 
integral ( pq) for calculating Q) ' should be replaced by the Debye-type finite radius D 
of screening for such system. Then we obtain: 

,3D 

The stress tensor and heat flux vector for obtained distribution function f\(xi;t) in 
the first approximation are: 

P (n; t)=P (n; t) 7 -ae(w (n; t))-2 V S (n; t) , (32) 

where P is the pressure, S ( r i',t) is the shift tensor of velocities, as is the coefficient of 
bulk viscosity 

as = -a 4 n 2 g2 (o + \n,j3^J V irrnkT, (33) 
7] is the coefficient of shear viscosity 

11 = I* + ^ T ( X + ^™^ 2 l3 )) B (0) • (34) 

The heat flux vector has the form 

q(r 1 ;t) = -X(V-T(r 1 ;t)), (35) 

where A is the thermal conductivity: 



3k 5 , 2kT ( 2 , / , . x\ t , , 
2m 4 V rn \ 5 v y / 

We can consider some particular expressions for the quantities as, rj and A. If n — > 0, then 
as, ?7, and A coincide with respective quantities obtained by solving Boltzmann's equation 
for low density systems of point-like charged particles If Z = 0, then we obtain the 
results of RET-theory [TL4fl . 



4 Numerical calculations 



The numerical calculations were carried out for the viscosities as (p3|) , rj ( p4|) and thermal 



conductivity A (|36|), where the dense once-ionized argon was chosen as a dense one- 
component plasma in a homogeneous neutralizing continuum. In the case Z = 0, the 
obtained results were compared with [2vJ. In these papers dense neutral argon was con- 
sidered. Its atoms were modelled as hard spheres, and obtained results are in a good 
agreement with |23| if the theory parameter a is chosen correctly. 



The binary correlation distribution function g2(a + \n, (5) was taken from [24j where it 
is presented as the functional of a and n, but not of (3: 



92 



(a, nj 



71 

1 no" 

12 



x 



n 

1 no" 

6 



(37) 



S 



can give for high densities incorrect values being less than a. Therefore it is presented in 



the form, that was proposed in |pq| , where the hard spheres diameter is also taken into 
account: 

A 



D 



1 



a- 



2A 



1 + AZea 



1 



- A y nm 

~2AJ VfcT 



1 



-i 



(38) 



where 



A = -rcna 3 . 
6 



Such substitution of a, D and (72(0"; n) allows to obtain transport coefficients as functions 
of density n and temperature T. The behaviours of se, rj and A were investigated for 



various values of n and T including regions mentioned in ||23|| , but only for the case giving 
value of D greater than a. The purpose of these calculations was to study the dependence 
of ae, 77 and A on the long-range potential when the density is high. As was expected the 
behaviour of such transport coefficients in these regions appears to be rather smooth and 
monotonous. But the small deviation is observed between our results and comparative 
ones when the long interaction is "switched on". Transport coefficients appear to be 
slightly sensitive to including the long-range potential. They decrease slowly being of the 
same order. These behaviour coincide with expectations. 

We note that the numerical calculation for ae, rj and A was carried out and compared 
with data from |^6[ ^7|. ^8[ . As one can see from the figures below, the good coincidence of 
both data is observed in the "high" -temperature region. Theory parameter a for different 
densities was borrowed from 



2S, BO, 31 



|. It is possible to improve the results by 
choosing a more precise but in general case a is the function of n, T, Z and finding for 
this dependence is a microscopic problem. 

As it can be seen from the figures below, the most essential deviations between of 
theoretical calculations from the experimental data arise in the low-temperature region 
for rj (figure [I]) and for A (figure 0). m re& l physical systems at low temperature one can 
approach to gas-liquid phase transition point. But our theory is not applicable within 
critical region near point of phase transition. Except that there is no phase transition in 
one-component system with only repulsive potential on the contrary to real systems. So 
our theory works well for description of transport processes at high temperatures far from 
phase transition point of corresponding real systems. 
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Figure 1: Temperature dependence of shear viscosity r] of neutral Ar at A = 0.1 (n = 
4.86 ■ 10 21 cm -3 ). Solid line represents results from theory, solid marked line represents 
data of |27j . Both r) and T are dimensionless. The transition relations to dimensional 
data read: r] dim = r\ ■ 10 -5 Pa ■ sec, T dim = T ■ 10 2 K. 
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Figure 2: Temperature dependence of thermal conductivity A of neutral Ar at A = 0.075 
(n = 3.644 ■ 10 21 cm -3 ). The legend is the same as for figure |l|. Both A and T are 
dimensionless. The transition relations to dimensional data read: A^ m = A - 10 -2 Wt/(m- 
A), T dim = T • 10 2 K. 
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Figure 3: Temperature dependence of thermal conductivity A of neutral Ar at A = 0.0125 



6.074 ■ 10 20 cm 



were taken from [p8 
figure 0. 



~ 3 ). The legend is the same as for figure [I], but experimental data 
The transition relations to dimensional data read the same as for 
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Figure 4: Temperature dependence of thermal conductivity A of once-ionized Ar at A 



0.0126 {n = 6.123 • 10 cm ). The legend is the same as for figure 
relations to dimensional data read: Xdim = A • 10 _1 Wt/(m ■ K), Tdi m - 



]. The transition 
T • 10 3 K. 
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